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Anisotropic shear melting and recrystallization 
of a two-dimensional complex (dusty) plasma 
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A two-dimensional plasma crystal was melted by suddenly applying localized shear stress. A stripe 
of particles in the crystal was pushed by the radiation pressure force of a laser beam. We found that 
the response of the plasma crystal to stress and the eventual shear melting depended strongly on the 
crystal's angular orientation relative to the laser beam. Shear stress and strain rate were measured, 
from which the spatially resolved shear viscosity was calculated. The latter was shown to have 
minima in the regions with high velocity shear, thus demonstrating shear thinning. Shear-induced 
reordering was observed in the steady-state flow, where particles formed strings aligned in the flow 
direction. 

PACS numbers: 52.27.Lw, 52.27.Gr, 82.70.Dd 
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I. INTRODUCTION 

Shear melting is a non-equilibrium process where a 
crystalline lattice is ruptured and eventually melted by 
an external shearing force. Shear melting and related 
phenomena were experimentally studied in soft con- 
densed matter systems such as colloids [1-3] and complex 
(dusty) plasmas [3|-l6|. 

A complex, or dusty plasma is a suspension of small 
solid particles in an ionized gas 0]. Particles get 
charged by collecting electrons and ions from the am- 
bient plasma and interact with each other via a screened 
Coulomb pair potential. When these interparticle forces 
are balanced by other forces (gravity, neutral and ion 
drag, ambipolar electric fields), the particles can form 
regular structures called plasma crystals. In the presence 
of gravity, a single-layer (two-dimensional, 2D) plasma 
crystal can form [4-6]. In a 2D plasma crystal, the in- 
plane interparticle pair potential can be well approxi- 
mated by the Yukawa potential. The shear modulus of 
such crystals is on the scale of 10 -13 N/mm [6| , and there- 
fore they can be easily manipulated, e.g. by the radiation 
pressure force of a focused laser beam [H, [§| . 

Using plasma crystals as model systems, shear melting 
was studied experimentally in Refs. 0-0 • Proliferation of 
supersonic (for shear sound) dislocations was suggested 
as a possible mechanism of shear melting [(| . Shear flows 
in melted complex plasmas and their transport proper- 
ties, e.g shear viscosity were studied in Refs. |4i fiol. UA\ . 
One important aspect of shear melting that has not been 
studied so far is the role of relative angular orientation 
of the crystal and applied shearing force. 

In this work, we studied the response of a 2D complex 
plasma crystal to a suddenly applied shear stress. The 
crystal was pre-oriented to align the closely packed rows 
of particles at 0° or 30° (two principal orientations) rela- 
tive to the shearing force direction. The force was applied 
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FIG. 1: (color online) Shear stress is applied to a two- 
dimensional plasma crystal by the radiation pressure force of 
a rastered laser beam. The laser force is applied in the stripe 
—2 mm < x < 2 mm in the positive y direction. Particle tra- 
jectories in the camera field of view are shown during 1.67 s 
of steady-state shear flow (corresponds to Piaser = 2.4 W in 
Figs. 1214]) . Analysis was performed in the region indicated by 
the dashed-line rectangle, S(k) was calculated in the smaller 
area indicated by the dotted-line square. Before the laser was 
switched on, the plasma crystal was rotated so that the closely 
packed rows of particles formed an angle = 30° (see inset) 
or 0° with the laser beam. 



by the radiation pressure of a focused laser beam. 



II. EXPERIMENTAL METHOD 
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Our experimental setup was a modified GEC (Gaseous 
Electronics Conference) rf reference cell [6[ . Plasma was 
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produced using a capacitively coupled rf discharge in ar- 
gon at 0.66 Pa. A single layer of dust particles was sus- 
pended in the plasma sheath of the lower rf electrode. 
The microspheres made of melamine formaldehyde had a 
diameter of 9.19 ±0.09 /im, a mass m = 6.15 x 10 -13 kg, 
and acquired an electric charge of Q = — 16 000 ± 1600e 
[l2j . The suspension included around 9000 particles and 
had a diameter of « 60 mm. The mean interparticle dis- 
tance in the center was A = 0.51 mm [measured from 
the first peak of the pair correlation function g(r)]. The 
Wigner-Seitz radius was a = (7m) -1 / 2 = 0.27 mm, where 
n is the areal number density of particles. This cor- 
responds to the screening parameter n = a/Xo = 0.5, 
where is the screening length [12j. The neutral gas 
damping rate was v = 0.77 s _1 . 

In our experimental conditions, the particle suspen- 
sion self-organized in a highly ordered triangular lattice, 
as indicated by the high value of the (local) coupling pa- 
rameter r = Q 2 (47re ak B T)- 1 ~ 9000, where T is the 
particle kinetic temperature. When this plasma crystal 
settled, we rotated it to align one of the closely packed 
rows at = 30° or 0° with respect to the laser beam, 
see inset in Fig. [TJ To rotate the crystal, we temporar- 
ily placed a small permanent magnet on the upper glass 
window [l3|| . After the crystal was properly aligned, the 
magnet was removed. 

Shear stress was applied to the plasma crystal by push- 
ing a stripe of particles by the radiation pressure force of 
a rastered laser beam, as shown in Fig. [TJ A Millen- 
nia PRO 15s J solid-state 532 nm laser with a maximum 
CW output power of 15 W was used. The laser beam 
coming at a grazing angle of ~ 8° was focused on the 
particle suspension to a spot smaller than the interpar- 
ticle distance and then rapidly scanned to draw a rect- 
angular stripe on the suspension. The scanning was per- 
formed by galvanometer-driven mirrors oscillating at a 
frequency of around 300 Hz. The particles within the illu- 
minated stripe reacted to the average radiation pressure 
force and were pushed in the direction of the laser beam. 
Shear stress was thus created in the particle suspension, 
its magnitude was controlled by varying the laser output 
power. 

The lattice layer was illuminated by a 660 nm, 20 mW 
diode laser, with its output expanded into a horizontal 
sheet by a cylindrical lens. A video camera (Photron 
FASTCAM 1024 PCI) was mounted above the chamber, 
capturing a top view of the particle suspension with a size 
of 42.7 x 42.7 mm 2 and resolution 1024 x 1024 pixel 2 . 
The camera lens was equipped with a narrow-band in- 
terference filter to admit only the illumination laser light 
scattered off the particles. The recording rate was set at 
60 frames per second. 



III. RESULTS 

We applied increasing levels of shear stress to the crys- 
tal. At low stress, the crystal deformed elastically. At 
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FIG. 2: Averaged particle velocity V y in the laser- illuminated 
stripe as a function of time, for two principal orientations of 
the plasma crystal. Different curves are for the manipulation 
laser power Pi aS er ranging from 0.2 W to 2.2 W, in increments 
of 0.2 W. The laser was switched on at t = and off at 
t = 10 s. 



higher stress, the following stages were observed: defect 
generation while in a solid state, onset of plastic flow, 
and fully developed shear flow. In this paper, we will 
study in detail the last two stages. 



A. Development of shear-induced melting 

For an overview of the succession of stages leading to 
shear melting, we measured the average velocity V y [1J] 
of all particles in the laser-illuminated stripe as a func- 
tion of time, see Fig. [2j For 6 = 0° (upper panel), there 
is a rather gradual change in the V y (t) curves for increas- 
ing laser power. The particle suspension becomes liquid 
and flowing at higher laser powers, but there is no evi- 
dent threshold for shear melting. In contrast, for = 30° 
(lower panel) the shear melting has a pronounced thresh- 
old character: By increasing the laser output power from 
1.8 W to 2.0 W, the V y (t) curve changes dramatically. 
For laser powers Pi aS er > 2.0 W, the applied stress ap- 
parently exceeds the yield stress and the lattice eventu- 
ally shear- melts. For Pi aS er = 1-6 — 1.8 W, the lattice 
defects (mostly dislocation pairs) are continuously gen- 
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FIG. 3: (color online) (a),(c) Particle trajectories during 
0.23 s and (b),(d) the lattice vorticity (V x v) z , where v 
is the velocity of individual particles, at the onset of plastic 
deformation. The particles in the shaded area were pushed 
by the laser (Pi aS er = 2.4 W), as indicated by the arrows. The 
lattice orientation before the laser was applied was (a),(b) 
= 0° and (c),(d) = 30°. The lattice vorticity reveals shear 
motion and is used here to visualize moving dislocations [TEj]. 



erated which however do not lead to shear melting on the 
time scale of the present experiment. The finite values of 
V v in this case reveal a kind of a creep flow in the lattice 

A 

Shear melting occurs more easily when closely packed 
rows of the crystal are aligned with the laser beam (0 = 
0°). This configuration favors dislocation nucleation and 
subsequent motion, since they glide in the direction of the 
laser beam 0|, see Figs.[3]^a),(b). For = 30°, the plastic 
deformation of the lattice occurs via dislocations gliding 
at a certain angle with respect to the laser beam, see 
Figs. [3(c) , (d) . Therefore, the melting happens at higher 
Pi aS er and has a more pronounced threshold character. 

The first peak at t ~ 0.5 s in Fig. [2] corresponds to 
in-cage oscillations of particles. Its presence is a mani- 
festation of a gradual crossover from elastic to plastic de- 
formation. It is naturally more pronounced for = 30°. 

We measured the displacement of the laser-illuminated 
particle stripe by integrating the average velocity V y [14[ 
of all particles in the stripe, see Fig. HJ After the laser 
is switched on at t = the displacement scales as ex t 2 , 
which indicates that the particles accelerate due to the 
action of the applied force. After c± 0.5 s, the particle 
motion slows down. As long as Pi aS er is below a threshold 
the curves tend to saturate, whereas above the threshold 
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FIG. 4: Displacement of the laser-illuminated particle stripe 
as a function of time, for two principal orientations of the 
triangular lattice. Different curves are for the manipulation 
laser power Pi aS er ranging from 0.4 W to 2.2 W, in increments 
of 0.2 W. The laser was switched on at t = and off at 
t = 10 s. 



the lattice shear-melts. Then y oc where a = 0.9 — 1.7, 
depending on the laser power. Notice the final push- 
back, even of a liquid complex plasma after the laser was 
switched off at t = 10 s. This is caused by the elastic 
response of unmelted lattice and by the pressure gradient 
in the liquid state. 



B. Shear-induced reordering 

We observed shear-induced reordering in a flowing liq- 
uid complex plasma. We made two related observations. 
First, a flowing liquid complex plasma is not completely 
isotropic. Rather, on top of isotropic particle distribu- 
tion appears a weak triangular order with = 0° even 
in the cases when the original triangular lattice was ori- 
ented at 30° (see Fig.[5j third line). Similar anisotropy in 
a liquid complex plasma was observed in Ref . [17j , where 
the initial triangular order with = 0° was not com- 
pletely destroyed in a shear melted state. Second, after 
the laser was switched off, the liquid complex plasma re- 
crystallized in a triangular lattice with = 0° (see Fig.[5j 
last line), even though the lattice was oriented at 6 = 30° 
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FIG. 5: Shear-induced reordering in a 2D complex plasma. 
The left column shows snapshots of particle positions, the 
right column shows corresponding static structure factor S(k) 
(all from the area indicated by the dotted-line square in 
Fig. [T}. Undisturbed triangular lattice was initially oriented 
at = 30°, first row. The lattice was distorted elastically 
t = 0.67 s after the laser was applied (Piaser = 2.4 W), sec- 
ond row. During the subsequent steady-state shear flow, the 
complex plasma was in a liquid state which is, however, not 
completely isotropic and has a nascent triangular order with 
= 0°, third row [S(k) averaged over t = 5 — 10 s]. After the 
laser was switched off, the lattice re-crystallized with — 0°, 
last row [S(k) averaged over t = 10 — 13.3 s]. 

before melting. 

A flowing liquid complex plasmas, therefore, favors 
particle ordering where "closely packed lines" are ori- 
ented along the flow. It was further observed that even 
outside of the shear flow the adjacent crystalline regions 
gradually change their orientation from = 30° to 0° 
and stay in that orientation after cessation of the flow. 

The observed particle ordering is similar to the for- 



mation of particle strings aligned in the flow direction in 
sheared 2D colloidal dispersions [18| and simulated liquid 
of soft disks [l9j. Here, we observed it for the first time in 
a sheared complex plasma. The underlying mechanism 
might be similar to the phenomenon of lane formation ob- 
served in binary systems where particles of one sort were 
driven through a cloud of particles of another sort, e.g. 
in colloidal dispersions [3|, [20|, [21( and complex plasmas 

asm. 

C. Rheology of liquid complex plasma 

In this section, we discuss the rheology of a steady- 
state shear flow in liquid complex plasma. The time t w 
for the flow to reach steady state varied depending on 
the initial crystal orientation and the laser power ap- 
plied, see Fig. [2j For example, t w ~ 4 s for 6 = 30° and 
flaser = 2.4 W. For this case, the particle trajectories in 
a steady-state shear flow are shown in Fig. [TJ The laser 
force is applied in the stripe —2 mm < x < 2 mm. Par- 
ticles far from the flow are only slightly disturbed and 
the original crystalline structure here is intact. There is 
a boundary layer between the flow and unmelted crys- 
tal. Profiles of time-averaged components of the par- 
ticle velocity v XyV (x) are shown in Fig. [6fa), and the 
longitudinal shear strain rate 7 = dv y /dx is plotted in 
Fig. [6^c) (solid curve). The profiles of the coupling pa- 
rameter are shown in Fig. [6fb), where r x%y corresponds 
to the particle kinetic temperature T X) y [24|. The values 
of r x ~ F y « 9000 far from the flow indicate a highly or- 
dered crystal. There is a small anisotropy in the coupling 
parameter, in particular at the edges of the flow, where 
T y < r xi indicating noticeable deviation of the system 
(driven in the y direction) from equilibrium. 

We note that the values of coupling parameter T y c± 
r x « 430 in the middle of the flow and F y « 300 at 
its edges are still above the equilibrium melting line in 
the (/€, r) coordinates as reported in simulations [25ll26|. 
This is because melting in our experiment is essen- 
tially a non-equilibrium phenomenon driven by external 
forces. Indeed, by comparing the relative flow velocity 
^fl ow ~ 7 A with the thermal velocity at the flow edges, 
vr y — y/kBTy/m, we conclude that these are about the 
same. Therefore, the effective "non-equilibrium" cou- 
pling parameter (in the y direction), which is inversely 
proportional to v| ow + v\ , is substantially reduced with 
respect to the "equilibrium" value of F y . 

To characterize the flow, we calculate the shear stress 
a xy in the particle suspension using the positions of indi- 
vidual particles, as explained below. We used a method 
where the stress at the position of a particle is built up 
from individual contributions of its neighbors. In this 
approach, the stress tensor is given by Q: 

1 2 f r®rdV 
where g(r) is pair correlation function and V(r) is the 
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FIG. 6: (a) Averaged particle velocity v x , y , (b) coupling 
parameter r x , y , (c) shear stress a xy and shear rate 7 = 
dvy/dx, and (d) kinematic viscosity v — p~ 1 a xy /j (where 
p — mn is the areal mass density) as functions of the trans- 
verse coordinate in the steady-state shear flow. The di- 
mensionless viscosity 77* = v / (cj p d a 2 ) is also shown, where 
w p d = [Q 2 /2neoma 3 ] 1 / 2 is the relevant plasma frequency and 
a — (7m) -1 / 2 is the Wigner-Seitz radius. The averaging 
was performed in narrow rectangular bins with the width 
of Ax = 0.42 mm during 5.7 s. The dotted line in (d) 
is a fit of experimental data with v — z/o(l — cry 2 ), where 
i/q — 5.92 mm 2 /s and a — 0.14 s 2 . The lattice was oriented 
at — 30° before the laser was switched on, the laser power 

WaS Pl aS er = 2.4 W. 



pair interaction potential. Then for our 2D system we 
obtain in polar coordinates: 

&xy = —-n 2 / drr 2 —— / d(j) cos cj) sin (j) g(r, (p). 
2 Jo dr Jo 

(2) 

To calculate a xy from Eq. (j2]), V(r) was approximated by 
the Yukawa potential and 16 shells of neighbors were in- 



FIG. 7: Pair correlation function g(r, eft) of complex plasma (a) 
in a crystalline area far from the shear flow, x = 8.1 mm, (b) in 
the middle of the shear flow, (c),(d) at its edges, x = —1.9 mm 
and 1.9 mm, respectively (see Fig. [6}. The peak amplitudes 
are up to 41 in (a) and therefore appear saturated in this 
color coding. Elliptically distorted g(r,</>) in (c),(d) indicates 
strong shear stress. Only 4 shells of neighbors are shown here, 
shear stress was calculated using 16 shells. 



eluded in g(r, (p). The result is shown in Fig.[6fc) (dotted 
curve). Note that we calculate only the potential part of 
shear stress, since the remaining kinetic part is negligible 
in our experimental conditions (T > 300) p7j . 

It is clear from Eq. (j2j) that the shear stress is de- 
termined by the sin 20 component of g(r,(fi). As it was 
shown in Ref. [17], the interparticle spacing A(c J >) [repre- 
sentative of the first peak of g(r, (/))] in a sheared liquid 
complex plasma indeed has a strong sin 2<p component. In 
Fig. [7J we present the complete g(r, 0), where the sin 20 
component shows up as an elliptical distortion particu- 
larly clearly seen in the first peak, see panels (c) and 

(d). 

Given the shear stress and strain rate, the kinematic 
viscosity of liquid complex plasma was calculated as 
v = p~ 1 (j xy / A f^ where p = mn is the areal mass den- 
sity. The result for v is shown in Fig. [6f d) (solid curve). 
Shear viscosity has two local minima at the edges of the 
flow, x = —1.9 mm and 1.5 mm, where the shear rate 
is maximal. This indicates shear thinning, i.e., a kind of 
non-Newtonian behavior of a liquid whose viscosity di- 
minishes for higher strain rates. Shear thinning in com- 
plex plasmas was demonstrated in experiments [lOj and 
simulations [28|. In the region —1.4 mm < x < 0.6 mm, 
where 7^0 the shear viscosity cannot be calculated as 
Vxy/i and therefore is not shown in Fig. [6fd). 

The found values of kinematic viscosity compare well 
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FIG. 8: Kinematic viscosity v is shown as a function of 7 2 , 
revealing shear thinning. The line is a fit of experimental data 
with v — ZA)(1 — cry 2 ), where z/o — 5.92 mm 2 /s and a — 0.14 s 2 . 
The dimensionless viscosity of — v/(u p da 2 ) is also shown. 

with the results of Ref. [4|. In the present work, v ~ 
2.5 mm 2 /s at the edges of shear flow where F = 300—370. 
In Ref. [J, v ~ 2 mm 2 /s, for similar values of F and 
k. In both experiments, shear thinning was apparently 
well pronounced. In Ref. [4], viscosity was calculated 
from the fits of experimental particle velocity profiles to 
the solution of a modified Navier- Stokes equation taking 
neutral gas drag into account. It is remarkable that the 
two very different methods of calculating shear viscosity 
give similar results. 

To characterize shear thinning, we plot kinematic vis- 
cosity v as a function of 7 2 , see Fig. [U Viscosity de- 
clines roughly three times when 7 2 increases from to 
5 s -2 . Note that simultaneously the coupling parameter 
F varies by 14 — 42%, see Fig. [6fb), due to the viscous 
heat release [29| and higher temperature at the edges 
of the flow, where the shear rate is higher. The varia- 
tions in 7, F , and v are self-consistent due to the depen- 
dence of shear viscosity on temperature [To|- The value 
of viscosity in the limit of low shear, = 5.92 mm 2 /s 
can be c omp ared to the results of equilibrium simula- 
tions [23, [30j . It corresponds to dimensionless viscosity 



77* v/(uj pd a 2 ) = 0.8, where uj pd = [Q 2 /27re raa 3 ] 1/2 
and a = {itn) -1 ' 2 . In the equilibrium 2D MD simula- 
tion of Ref. [27j , where viscosity was calculated using the 
Green-Kubo formula, 77* = 0.5 — 1.6 for the highest re- 
ported values of r = 125 — 130, for a comparable value 
of K. 



IV. SUMMARY 

Shear induced melting of a 2D plasma crystal was stud- 
ied experimentally. Shear stress was applied by a focused 
laser beam rastered over a rectangular area in the crystal. 
For the first time, the study was angle-resolved, i.e., the 
closely packed rows of the crystal were oriented at 6 — 0° 
or 30° (two principal orientations) with respect to the 
shearing force. Shear melting was found to be strongly 
anisotropic. The transition of the crystal to flowing liq- 
uid with increasing shear stress was rather gradual for 
6 = 0° and had a distinctly threshold character for 30°. 
This is explained by the role of dislocations in the crys- 
tal plastic deformation. Their nucleation and motion was 
facilitated for 0° and retarded for 30°. 

The rheology of the steady-state shear flow in liquid 
complex plasma was studied. Shear stress and strain rate 
were measured with spatial resolution and from these, 
the shear viscosity was calculated. Shear viscosity had 
minima in the regions with high velocity shear, demon- 
strating shear thinning. 

Shear-induced reordering was observed in flowing liq- 
uid complex plasma, where particles tended to form 
strings aligned in the direction of the flow. Different 
from laning in driven binary complex plasmas, forma- 
tion of strings in sheared complex plasma was observed 
for the first time. 
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